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Q ■ Abstract 



. The first, to our knowledge, calculation of neutrinoless double beta decay 

£3 ■ (Oz//3/3-decay ) matrix elements within the self-consistent renormalized Quasi- 



X 



particle Random Phase Approximation (SRQRPA) is presented. The con- 
tribution from the momentum-dependent induced nucleon currents to 
decay amplitude is taken into account. A detailed nuclear structure study 
includes the discussion of the sensitivity of the obtained SRQRPA results for 
0^/?/3-decay of 76 Ge to the parameters of nuclear Hamiltonian, two-nucleon 
short-range correlations and the truncation of the model space. A compari- 
son with the standard and renormalized QRPA results is presented. We have 
found a considerable reduction of the SRQRPA nuclear matrix elements, re- 
sulting in less stringent limits for the effective neutrino mass. PACS numbers: 
23.40.Hc, 21.60.Jz, 27.50.+e, 27.60.+j 
The neutrino mass and mixing of neutrinos are the main subject of the elementary 
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particle physics nowadays. The experiments looking for oscillation of solar, atmospheric and 
terrestrial (LSND-experiment) neutrinos constitute evidence for a new physics beyond the 
Standard Model JTJ. The current constraints imposed by the results of neutrino oscillation 
experiments allow to construct the spectrum of mass of neutrinos |2[. The predictions differ 
from each other by different input, structure of neutrino mixing and assumption, e.g., on 
the phases and fundamental character of neutrinos. 

An important quantity to limit the space of possible neutrino mixing schemes is the 
effective Majorana neutrino mass 



(mee) 



(U ck ) 2 m k r]£ 

k 



CP 



(1) 



where U e k, vrik and r/J? p are unitary mixing matrix elements, mass eigenstates and the relative 
CP-phases of neutrinos, respectively. The value of (m ee ) depends on the specific model of 
neutrino mixing and its predictions fall in the range of 10~ 3 eV to few eV ||. The upper 
bound on (m ee ) can be deduced from the current experimental lower bound on the half-life 
of neutrinoless double beta decay (Oi//?/?-decay) as follows 

(mee) < m \ (2) 

l M (m oc }IV J l/2 Go1 



Here, m e , Goi and Mi m \ are the mass of electron, the integrated kinematical factor PlJlC 
and the Oz/^/S-decay nuclear matrix element, respectively. 

The present most stringent lower bound on T^ -63 ^ has been measured for 76 Ge by 
Heidelberg-Moscow group and is equal to 5.7 x 10 25 years |||7|. A condition for obtaining 
reliable limit for fundamental particle physics quantity (m cc ) is that the nuclear matrix ele- 
ments governing the Majorana mass mechanism of 0z//3/3-decay can be calculated correctly 
0. However, the practical nuclear structure calculation always involves some approxima- 
tions, which make it difficult to obtain an unambiguous limit on (m ee ). We note that the 
difference between the previous (M^ mc ^ = 4.18, i.e., (m ee ) < 0.18) [|J and the recent more 
advanced (M, m \ = 1.67, i.e., (m ee ) < 0.46) [§] shell model calculations is significant. There 
is also an other group of nuclear structure calculations which include the proton-neutron 



Quasiparticle Random Phase Approximation (pn-QRPA) || and its extensions Some 
of them suggest the upper bound on (m ee ) to be in the 0.1 eV range ||. 

The aim of the present letter is to discuss the nuclear physics aspects of the pn-QRPA 
0, the renormalized pn-QRPA (pn-RQRPA) [L2|,[13| and the self-consistent pn-RQRPA 
(pn-SRQRPA) p3| , [TT| calculation of the nuclear matrix element M, m \ for 0i//3/3-decay of 
76 Ge. We note that the pn-SRQRPA results on Oz//3/3-decay matrix elements have been not 
presented in the literature till now. 

We shortly present the main differences between the above mentioned three QRPA ap- 
proaches. 

The pn-QRPA has been the most popular theoretical tool in description of the (3 and 
(3(3 decays of medium and heavy open shell nuclei. However, the pn-QRPA develops a 
collapse beyond some critical value of the particle-particle interaction strength of nuclear 
Hamiltonian close to its realistic value. This phenomena makes predictive power of the 
obtained results questionable. 

By implementing the Pauli exclusion principle (PEP) in an approximate way in the pn- 
QRPA one gets the pn-RQRPA [[TJ,|rjJ, which avoids collapse within a physical range of 
particle-particle force and offers more stable solution. This fact has been confirmed also 
within the schematic models. In addition it was found that by restoring the PEP a better 
agreement with the exact solution is obtained ||14j| . 

The selfconsistent pn-RQRPA (pn-SRQRPA) |T5|JT7| , p!6[| , which is more complex version 
of the RQRPA, is becoming increasingly popular to describe strongly correlated Fermion 
systems. The pn-SRQRPA goes a step further beyond the pn-RQRPA. In the pn-SRQRPA 
at the same time the mean field is changed by minimizing the energy and fixing the number 
of particles in the correlated ground state instead of the uncorrelated BCS one as it is done in 



the other versions of the QRPA (pn-QRPA, pn-RQRPA) |IJ],0,[T6| . Thus the pn-SRQRPA 
is closer to a fully variational theory. 

We proceed by writing the expression for the Oz//3/3-decay nuclear matrix element. We 
note that the contribution to M^ mcc ^ from the induced pseudoscalar term of the nucleon 



current were not considered before. Recently, it has been found that it is significant and 



leads to a modification of Gamow- Teller and new tensor contributions of M^ mc ^ ||11|| . Thus 



M( mcc ) is given as sum of Fermi, Gamow- Teller and tensor contributions 

M F 

M {mcc) = f + M GT + M T (3) 

9a 

with 9a = 1.25. Expressed in relative coordinates and using the second quantization for- 



malism M{ mcc ) takes the form 



jp jn J 
jn' jp' J 



x(p(l),p\2);J\f(r 12 )r^O {mcc) (l2)f(r 12 )\n(l),n\2);J) 

x(0+ || [cpf n ,]j || rm f )(J*m f \rmi){rmi || [c+c„]j || 0+). (4) 

Here, f{rvz) is the short-range correlation function and 0( mec )(12) represents the coordinate 
and spin dependent part of the two-body Oz//3/3-decay transition operator 

<mcc) (12) = -M^l + H GT (r 12 )a 12 + H T {r 12 )S 12 , (5) 
9a 

where a\ 2 = <7i ■ a 2 , S 12 = 3(<?i • qa 2 ■ q) - a\ 2l and 

H K (r 12 ) = A" T -^L fe,^) dq, (K = F, GT, T) (6) 
T^9A r vi^ q + E m (J) — (E l + Ej )/2 

with 

h F {q 2 ) = gW)9A, 

h (,?\ n 2/ 2x | l ffp(g 2 )g 4 2 g A (q 2 )gp(q 2 )q 2 
h GT (q ) - 9a(Q ) + 3^2- " 3 ^ ' 

kT{q ] = 3 2^ 3^T" (?) 



Here, i? is the mean nuclear radius | IT| . E l , E* and E m (J) are respectively the energies of 
the initial, final and intermediate nuclear state with angular momentum J. The momentum 
dependence of the vector, axial- vector and pseudoscalar formfactors (gv(q 2 ), fiu(<? 2 ) and 



9p(q )) is given in Ref. JTT 



We note that the first sum on the r.h.s of Eq. ([|) represents summation over all multipo- 
larities. The form of the one-body transition densities [see Eq. (Q)] to excited states \J n rrii) 
and | J n mf) generated from the initial (A,Z) and the final (A,Z+2) ground states \0f) and 
Ot) within the pn-RQRPA and the pn-SRQRPA is the same and can be found together with 
other details of the nuclear structure in Refs. [ P~5] , [TT] , |TT] | . The difference consists only in the 
calculated value of the renormalized coefficients D, which in the case of the pn-QRPA is 
just equal to unity. The overlap factor entering the expression (f|) can be find in Ref. [|18[ . 

The calculation of 0z//3/3-decay matrix elements of 76 Ge are performed within two model 
spaces both for protons and neutrons as follows: i) The model space I (m.s. I) consists of 
the full 3 — 4huj major oscillator shells and has been considered in the pn-QRPA studies of 
Ref. (9- levels model space), ii) The model space II (m.s. II) comprises the full 2 — 5hu 
major shells (12 levels model space). 

The single particle energies were obtained by using a Coulomb-corrected Woods-Saxon 
potential with Bertsch parameterization. Two-body G-matrix elements were calculated 
from the Bonn one-boson exchange potential within within the Brueckner theory. In the 
pn-QRPA and the pn-RQRPA approaches pairing interactions have been adjusted to fit the 
empirical pairing gaps. In the pn-SRQRPA approach the pairing matrix elements of the NN 
interaction have not been rescaled as the mean field is directly related to the excited states. 
The particle-particle and particle-hole channels of the G-matrix interaction of the nuclear 
Hamiltonian H are renormalized by introducing the parameters g pp and g p ^, respectively. 

In Fig. |l] the calculated partial matrix elements Mp, Mqt and Mp are plotted as 
function of g pp (g p ^ = 1.0) for the pn-RQRPA and pn-SRQRPA approaches. The larger 
12-levels model space is considered. One find a strong dependence of Mp and Mqt on g pp . 
The smallest Mp matrix element is rather insensitive to this parameter. In general, the 
behavior of plotted matrix elements is similar for both approaches. Nevertheless, the Fermi 
and GT matrix elements of the pn-SRQRPA one reach zero value inside the physical range 
of g pp parameter (0.8 < g pp < 1.2). 

In Fig. |2| the multipole decomposition (according to intermediate multipoles J n ) of 



matrix element M( m \ for the pn-RQRPA and pn-SRQRPA approaches and two values of 
9pp (dpp — 0.0 and g pp = 1.0) is presented. The model space and g p ^ is the same as in Fig. 
[1]. The filled and open bars represent the calculations with and without consideration of the 
two-nucleon short-range correlations (s.r.c), respectively. One clearly sees that the s.r.c. 
suppress the contributions from higher multipolarities. It just proves that the Majorana 
neutrino mass exchange mechanism is predominantly long-range interaction in the nucleus. 
We notice also that the contributions associated with even parity multipolarities are more 
influenced by the change of g pp (especially in the case of the pn-SRQRPA) as those associated 
with the negative ones. The differences between the even parity multipole contributions 
calculated via the pn-RQRPA and the pn-SQRPA approaches is remendeous for g pp = 1.0. 
This effect can be attributed to the different structure of the mean field in both approaches. 

In Fig. ^| the nuclear matrix elements M( mcc ) as function of particle-particle strength g pp 
for all three QRPA approaches (pn-QRPA, pn-RQRPA, pn-SRQRPA) is shown. One notices 
that the pn-RQRPA and the pn-SRQRPA results are much less sensitive to the truncation 
of model space and g pp as the pn-QRPA ones. The unpleasant feature of the pn-QRPA and 
the pn-SRQRPA is that for the larger model space M( mcc ) crosses zero within the physical 
interval of g pp what increase uncertainty in respect to the extracted value on (m ec ) from 
non-observability of 0i//?/3-decay considerably. 

The real value of the g pp parameter is expected to be close to unity. In Table [I] we 
present M( mcc ) for this value of g pp parameter together with the corresponding upper bound 
on (m ee ) deduced from the Heisenberg-Moscow experiment. One notices that only the 
pn-QRPA results are sensitive to the truncation of model space. The new calculated pn- 
SRQRPA results predict considerably larger value for upper limit on (m ee ) in comparison 
with the pn-QRPA and the pn-RQRPA ones. It is due to the mean-field modifications by 
the number self-consistent approach which use to diminish the pairing gap |L9| . 

In summary, the nuclear matrix element associated with the light Majorana neutrino 
mass mechanism of 0u(3f3-decaj of 76 Ge has been calculated in the selfconsistent pn-RQRPA 
approach which received much attention recently. This approach predicts considerably sup- 
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pressed Oz/^/S-decay matrix element, what, when confronted with the present experimental 
data, leads to a less stringent limit on the effective Majorana neutrino mass parameter 
(m ee ). The reduction of the matrix element is due to the modification of the mean field 
by the combined selfconsistent BCS+RQRPA approach which influences mostly the even 
parity multipole contributions. 

This work was supported in part by the State Committee for Scientific Researches 
(Poland), Contract No. 2 P03 B00516. 
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TABLES 

TABLE I. Calculated nuclear matrix element Mi m \ for O^/3/3-decay of 76 Ge within the 
pn-QRPA, pn-RQRPA and the pn-SRQRPA and the corresponding upper limit on the effective neu- 
trino mass (m ee ) deduced from the experimental lower bound for a half-life Ty^{ Ge) > 5.7 x 10 25 
years Q. Here, m.s. I and m.s. II denote the model spaces which comprice the full 3 — 4huj and 
2 — 4huj major oscillator shells, respectively. 



pn-QRPA pn-RQRPA pn-SRQRPA 



9 P h 


m.s. I 


m.s. II 


m.s. I 


m.s. II 


m.s. I 


m.s. II 


0.8 


2.76 


1.66 


3.15 


2.80 


0.72 


0.65 




0.28 


0.46 


0.24 


0.27 


1.1 


1.2 



M, 



(m ee ) 



1.0 



2.54 
0.30 



1.53 
0.50 



2.95 
0.26 



2.61 
0.29 



0.66 
1.2 



0.59 
1.3 
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FIGURES 
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FIG. 1. The Fermi (Mp, dashed line), the Gamow- Teller (Mgt, solid line) and tensor 
(Mt, dot-dased line) nuclear matrix elements for the 0u(5(3-decay of 76 Ge calculated within the 
pn-RQRPA (a) and pn-SRQRPA (b) are plotted as a function of the particle-particle coupling 
constant g pp . The 12-level model space II (the full 2 — 4hu major oscillator shells) was assumed 
and g ph = 1.0 
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FIG. 2. The contributions of different multipolarities to the nuclear matrix element M^ mee ^ 
for the Oi//3/3-decay of 76 Ge calculated within the pn-RQRPA (left figures) and pn-SRQRPA (right 
figures) without (open bar) and with (closed bar) consideration of the nucleon short-hand correla- 
tions (s.r.c). The model space is the same as in Fig. 1 and g p h = 1.0. The results in upper and 
lower graphs were obtained for g pp = 0.0 and g pp = 1.0, respectively. 
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FIG. 3. The nuclear matrix element M(, mca ^ for the O^/3/3-decay of Ge calculated within the 
pn-QRPA (dashed line), pn-RQRPA (dot-dashed line) and pn-SQRPA (solid line) in the framework 
of the 9-level model space (a) (the full 3 — 4hoo major oscillator shells), and the 12-level model space 
(b) (the full 2 — major oscillator shells) is plotted as a function of the particle-particle coupling 
constant g pp . g p h was taken to be unity. 
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